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Toroidal metamaterial 
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It is shown that a new type of metamaterial, a 3D-array of toroidal solenoids, displays a significant 
toroidal response that can be readily measured. This is in sharp contrast to materials that exist in 
nature, where the toroidal component is weak and hardly measurable. The existence of an optimal 
configuration, maximizing the interaction with an external electromagnetic field, is demonstrated. 
In addition, it is found that a characteristic feature of the magnetic toroidal response is its strong 
dependence on the background dielectric permittivity of the host material, which suggests possible 
applications. Negative refraction and backward waves exist in a composite toroidal metamaterial, 
consisting of an array of wires and an array of toroidal solenoids. 
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Toroidal moments were first considered by Ya. 
Zel'dovich in 1957 They are fundamental electro- 
magnetic excitations that cannot be represented in terms 
of the standard multipole expansion [21] . The properties 
of materials that possess toroidal moments and the clas- 
sification of their interactions with external electromag- 
netic fields have become a subject of growing interest 
I, i i i, 1 0, i, i 0, [O, m, III Q 111 • In particular 
the unusual properties of non-radiating configurations, 
based on toroidal and supertoroidal currents, have been 
discussed d, d, [13] ■ Non-reciprocal refraction associated 
with an effective Lorentz force acting on photons prop- 
agating in a toroidal domain wall material has recently 
been predicted [7|. Whereas the importance of toroidal 
moments in particle physics has been established some 
time ago (see and the references therein), there is 
no known observation of toroidal response in the clas- 
sical electromagnetism Q. This is because the effects 
associated with toroidal response in materials that exist 
in nature are weak 0, S, i]. At the same time mod- 
ern technology allows arrays of toroidal moments with 
sub-milimeter particle-size to be manifactured and char- 
acterized Ts*]. Therefore, a metamaterial specifically de- 
signed to maximize the toroidal component of an interac- 
tion with an external electromagnetic field may provide 
an outcome that can open up new possibilities for obser- 
vation of toroidal response and its applications. 

The aim of the present paper is to provide an esti- 
mate for the magnitude of the response of an array of 
toroidal moments to an external electromagnetic field. 
The model is valid in a quasistatic regime (the size of 
each metaparticle is much smaller than the wavelength) 
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and with the assumption that a uniform current is in- 
duced in each metaparticle. It is shown that a toroidal 
metamaterial - an array of toroidal solenoids - shows a 
significant electromagnetic response and has intriguing 
electromagnetic properties. A characteristic feature is 
the strong dependence of the response bandwidth on the 
dielectric properties of the host material, which suggests 
possible applications. Backward waves and negative re- 
fraction exist in a composite toroidal material consisting 
of an array of wires and an array of toroidal metaparti- 
cles. 

It should be emphasized that the term "anapole" (of- 
ten used in nuclear physics problems) and the term 
"toroidal moment" refer to physically different objects 
(see the relevant discussion in [l3|). Indeed, an anapole 
does not radiate electromagnetic energy and does not in- 
teract with external electromagnetic fields. In contrast, 
toroidal moments interact with time- dependent electro- 
magnetic fields and radiate (scatter) electromagnetic en- 
ergy. Hence, an array of toroidal moments can act as 
an electromagnetic medium in the same way in which an 
array of split-ring resonators does. 

The discussion begins by considering a metaparticle in 
the form a toroidal solenoid as shown in Fig. [1] A con- 
ducting wire is wrapped around the toroid so that there 
are N turns between the entrance and the exit points 
for the uniform current fiowing through the wire. The 
two ends of the wire are connected through a capacitor 
C and through an additional loop in the equatorial plane 
of the structure to compensate the z-component of the 
magnetic dipole moment associated with the azimuthal 
component of the current fiowing in the toroidal section 
of the winding. An alternative way to eliminate the lat- 
ter is to use double winding, as shown in Fig. [if 6) p^ . 
It will be shown that the magnetic quadrupole moments 
of both the structures shown in Fig. [T]are equal to zero. 
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FIG. 1: (a) Toroidal moment - a toroid is wound continuously 
with wire and an extra loop of radius Ri = + 0.5(R/d)^, 
compensating the magnetic dipole moment of the toroidal sec- 
tion is introduced in the equatorial plane. L is the inductance 
of the solenoid and 5R is its loss resistance. C is an externally 
introduced capacitor. The origin of the co-ordinate system 
is the geometric center of the toroid. Sk ~ S{Ark/\Ark\), 
where S — -kB? and |rfc| = d. (6) A double- wound toroidal 
solenoid has a zero net magnetic dipole moment, (c) d and 
R axe the larger and the smaller radii of the toroidal surface 
(a/x^ -f — d)^ + = R^, respectively. 



Indeed, the toroidal dipole moment of a current density 
distribution j{r,t) is a vector given by [T^ [T3| 



^ J [{j.r)r-2r^j] dV. 



(1) 



On the other hand, the magnetic quadrupole moment of 
the same system is the second-rank tensor [T^j flst 



Q = ^ / [irxj)r + r{rx j)] dV. 



(2) 



A toroidal coil, assumed to be made of infinitely thin wire 
can be generated using the following equations 



X — {d — R cos Nip) cos ip, 



y = {d — R cos Nip) sin ip, 



Rs'mNip. 



(3a) 



(3b) 



(3c) 



If r = {x,y, z) is the radius-vector of an arbitrary point 
along the winding then ip is the angle between the x-axis 
and the projection of r onto the xy-plane. The assump- 
tion of a uniform current / flowing along the wire leads 
to the current element jdV being defined as jdV — Idr. 



Hence, given that dr = ^dtp^ equations ^ and ([2]) be- 
come 
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dip, (5) 



respectively. Evaluating ^ from ([3]), performing the in- 
tegrations in ([4]) and ^ and assuming iV > 1 yields. 



and 



irNIdR^ 
T = n 



Q = o, 



(6) 



(7) 



where n = (0,0,1) 15]. As dS]), ^ clearly show, the 
quadrupole moment of a toroidal solenoid is zero whereas 
its toroidal moment is non-zero. Although, as ([T]) and 
^ suggest, T and Q are measured in the same units 
[^m'^], they are physically different and this is reflected 
in their properties. For example, a DC toroidal moment 
(/ — const{t)) does not generate magnetic field outside 
its own volume [T5| . whereas a DC quadrupole moment 
generates non-zero magnetic field everywhere. In addi- 
tion, the radiation pattern of an AC toroidal moment is 
the same as that of an electric dipole [13, UBl , which is 
not the case for a magnetic quadrupole moment. The 
magnetic dipole moment of the toroidal coil 
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and, hence, this magnetic moment can be compensated 
by a loop of radius i?i = d\J\ -\- 0.5{R/d)'^ in the equato- 
rial plane of the torus, carrying current uniform current 
— /, as shown in Fig. [IJa). Note, that the magnetic 
quadrupole moment of this loop is zero and, hence, the 
magnetic quadrupole moment of the coil shown in Fig. 
[Ija) is identically zero. 

A toroidal coil, wound in a direction opposite to that 
of the coil defined by the system ([3]) relies instead upon 
the equations 



x' = {d — R cos Nip) sin ip, 



'/ = {d — R cos Nip) cos ip, 



z' = RsinNip, 



(10a) 



(10b) 



(10c) 



which are obtained from ([3]) by formally exchanging x 
and y in the latter. The toroidal, magnetic dipole and 
magnetic quadrupole moments are now T' = T , m' = 
—m and Q' — 0, so that it can be concluded that a 
double-wound toroidal solenoid has a zero net magnetic 
dipole and magnetic quadrupole moments, whereas the 
toroidal moment of the structure is twice that of a single- 
wound torus (Fig. [Tja). The possibility of eliminating 
the magnetic dipole moment of a toroidal coil by either 
adding a loop in the equatorial plane of the torus, or using 
a double winding solenoid has been discussed earlier [l^ . 

If either of the metaparticles shown in Fig. [l]is placed 
in an external inhomogeneous magnetic field B{r) and 
the characteristic length scale of variation of the field is 
much larger than the size of the metaparticle, then the 
leading-order response to the external field comes from 
the toroidal moment. Indeed, since the current / is uni- 
form there is no charge accumulation and the electric 
multipole moments are zero. The magnetic dipole mo- 
ment of the system can be elliminated as shown in Fig. 
[T]and the magnetic quadrupole moment is zero, as it has 
been demonstrated. Therefore, under these conditions it 
is toroidal moment T that generates the leading-order 
electromagnetic response of the system 0, [l^ [ij]. The 
net magnetic flux ^ in the solenoid can be written in the 
approximate form 



N 



(11) 



fe=l 



where Bk = B{rk) is the flux density in the center of 
a given loop. Since Sk is a vector pointing along the 
normal to this loop, then Sk can be approximated to 
Sk « S{Ark/\Ark\), where S = 7ri?^ provided that N 
is sufficiently large. Now \Ark\ ~ 2TTd/N, hence the flux 
(1) is 



SN 



N 



fc=i 



SN 
27rd 



B.dr. 



(12) 



dS 




In ([12)) the discreet distribution of loops has been ap- 
proximated with a continuous one, which is valid for a 
sufficiently large N and dS is a circle of radius d cen- 
tered at the origin. Using Stokes' theorem (|12p can be 
rewritten as 



SNd 



{n.[curl B]r= 



(13) 



FIG. 2: (a) Toroidal metamaterial - array of toroidal solenoids 
embedded in an isotropic dielectric material with dielectric 
permittivity e. (b) Composite toroidal metamaterial - array 
of toroidal solenoids and an array of long wires, parallel to 
the axis of the toroids. 



where once again the assumption for a weak magnetic 
field variation has been used. In (|13p [curl B]r=o stands 
for curl B evaluated at r = 0. 

The net energy VF of a toroidal metaparticle in an ex- 
ternal magnetic field is the sum of the contributions of 
each loop, hence 



N 



w 



-I^^Bk-Sk, 



(14) 



k=l 



where / is the strength of the current flowing in the 
solenoid. From andjM]) W takes the stan- 

dard form W = -T.curlB 0, [HI, i, [ll|, where T is 
given by ([5]). Suppose that the external electromagnetic 
field acting on the metaparticle is that of a monochro- 
matic plane wave {E, H) oc exp(— iwt + ik.r), where E 
is the electric field, H is the magnetic field, uj is the an- 
gular frequency and k is the wavevector. Kirchhoff's law 
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FIG. 3: The unit cell of a toroidal metamaterial. The mini- 
mum volume required for each toroid is 8R(d + Rf . 



implies that 



.dl_ 
''dt 



Q 
C 



d-P' 
~dt' 



(15) 



In L and are the inductance and the resistance of 
the solenoid, respectively, Q is the electric charge accu- 
mulated by the capacitor and ^' is the flux p3p . gener- 
ated by the local (microscopic) field H' . The complex 
amplitude of the current / is, therefore, 



/ = 



ILUTjlQ 



ujL) 



[curl H').n, 



(16) 



where r = {SNd)/2 and fio is the vacuum permeabil- 
ity. Consider now a toroidal metamaterial, a 3D array 
of toroidal metaparticles (Fig. ^ embedded in a back- 
ground medium. Two different types of host media are 
considered: ordinary isotropic dielectric media and ar- 
rays of long wires. The former situation is referred to as 
a "toroidal metamaterial" and the latter - a "composite 
toroidal metamaterial". It should be emphasized that if 
the wires do not cross the toroids, which is the situation 
depicted in Fig. [D(b), no interaction between the neigh- 
boring wires and toroids is possible. This outcome flows 
directly from the expression for the interaction energy of 
a toroidal moment with an external field. Indeed, since 
W = —T.curlB = —iioT.j, where j is the source of the 
field B, it is clear, that W is zero unless the locations of 
the toroidal moment and the field source coincide. The 
external field induces in each solenoid a time- dependent 
toroidal moment T that can be obtained from ^ and 
(fT6)) . In a similar way to the description of an array of 
magnetic dipoles in terms of its effective magnetization, 
an array of toroids can be characterized by an effective 
toroidization vector & — kT where k is the number of 
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FIG. 4: Dependence of the toroidal frequency lot on the ge- 
ometry of the unit cell. The transverse size R + d is fixed and 
the ratio R/d is the free parameter. 



toroids per unit volume. The macroscopic magnetization 
M of the material is then [3, 0, H, |Tl| 



M = V X 0. 



(17) 



Using ([6]), p6|) and the Lorenz-Lorentz formula H' 
H + M/3 allows (fT7)) to be written as 



M 



(fc X H).n 



1 + 



OJn 



(k X n), (18) 
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where the effective "toroidal" frequency uJj 
Lc^ (Hqt^ k)^^ has been introduced, ujq — (LC)^^^^ is the 
resonant frequency, 7 = SR/i is the damping factor, c is 
the speed of light in vacuum and fc^ = /c^ — (fc.n)^. The 
toroidal frequency plays a role similar to that of the ef- 
fective "plasma" frequency Wp in metamaterials involving 
arrays of wires [l3| ■ Equation (fT8|) clearly shows that the 
metaparticles are not responding like magnetic dipoles. 
Indeed, the magnetic flux density vector in the consti- 
tutive relationships is not proportional to the magnetic 
field locally, as is the case with metamaterials involving 
split-ring resonators (SRRs) [l3,[li,[il. As ^ shows, 
in order to maximize the electromagnetic response of the 
material one needs to minimize the toroidal frequency 
lut- Given R and d the maximum value of the number 
density is k = {8R{d+ R)'^)^^ (see Fig. [3]) and this corre- 
sponds to each toroid in Fig. [2] being in contact with its 
nearest neighbors. The toroidal frequency then is given 

by 



2 16c2(d + i?)2 

UJrp — 



TT^Rd^ 



(19) 



and it does not depend on the number of windings N. 
Note, that in deriving (HI]) L = fioN'^R'^ /{2d) has been 
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used [Ts'l. If the transverse size of the unit cell 2{R + d) 
is fixed (e.g. by the manifacturing process), the ratio 
R/d remains a free parameter and it is easy to show that 
(jJt reaches a minimum ujtq = 16c/{7Td^/3) at d = 3R 
and this value represents a fundamental upper limit for 
the strength of the toroidal metamaterial response in the 
quasistatic regime, given the transverse size of the unit 
cell 2{R + d). As Fig. |4] shows the function wt/ojto 
depends weakly on d/R in a relatively broad range of 
d/R values. However, structures with R << d are dis- 
advatageous since ujt ^ oo as R 0. It should be 
noted that the inductance of the solenoid is given by 
L = iiqN"^ R^ / (2d) provided that two conditions are met: 
(i) the torus is closely wound, i.e. N is sufficiently large 
and (ii) R « d. The first condition has already been 
imposed in deriving (fT2|) and, hence, its use in (fTO]) does 
not impose any further restrictions. Excellent agree- 
ment with experimental measurements has been obtained 
at = 15 13]. The inductance of a closely wound 
toroidal solenoid with arbitrary R and d is given by 
L' = noN"^ {(P - - i?2) f2^. Note, however, that 
at d = 3i?, which is the physically interesting param- 
eter range, the relative difference between L and L' is 
\L' - L\/L < 0.03. Thus, the use of L' instead of L in 
(fT9)) results in shifting the minimum of ujt from d = 3R 
to d« 3.2i?. 

With d =5 mm and R —1.7 mm the value of the 
toroidal frequency is iotqI2'k =28.3 GHz. In contrast, the 
micro-structures reported in ^is'l have a;To/27r «380 GHz. 
As ([T^ suggests, the validity of the quasistatic approx- 
imation requires lo « lot- This is because LOx^^d/c w 3 
and for this reason in what follows lo < lot/W is used. 
In addition, numerical simulations show that the lead- 
ing order toroidal moment T, given by Eq.®, describes 
the electromagnetic properties of toroidal solenoids with 
kd « 0.2 with a high degree of accuracy 10]. The host 
media, in which the toroids are embedded is modeled 
with the effective permittivity tensor 



(^d + {ep — l)nn 



(20) 



where 5 is the unit matrix and Sp = 1 pertains to the 
case of the host material being an isotropic dielectric. 

The combination e = 1 and Sp — 1 — (^j/cV-k^ ' where kp 
is the effective "plasma" wavenumber [21[, corresponds 
to a wire medium used as a host material. Note that at 
kz = (wave propagation perpendicular to the wires) Ep 



reduces to 



1 p 



For an electromagnetic wave propagating in the bulk 
sample of the toroidal metamaterial sketched in Fig. [5] 
the electromagnetic field is that of an "extraordinary" 
wave, i.e. E = {Ea:,0,E:,), H = {0,Hy,0) and k = 
(fcj,,0, fcz). The dispersion equation in this case is 
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FIG. 5: Influence of the Lorentz-Lorenz local field correction 
on the dispersion properties of a wave propagating in a bulk 
toroidal metamaterial at fe_Ln and E\\n. The full curves are 
obtained with H' — H + M /3 and the dashed curves - as- 
suming H' = H (low-density approximation). Real "Re" and 
imaginary "Im" parts of the normalized wavenumber kc/uj. 
(a) Toroidal metamaterial with 7/cJo = 5.5 x 10~^. (b) Com- 
posite toroidal metamaterial with •j/loo — 5.5 x 10"'' and 
LOp/LOo = 1.02\/2. In both (a) and (b) ujo/ujt = 0.07 and 
LOO = 2GHz. 



where 



4 

LO esr, 



l^LU 



(22) 



CO 



kUl + (3) 



t-ph.^ 



(21) 



If the frequency lo is sufficiently far from the resonant fre- 
quency loq the contribution of the toroidal component can 
be neglected (i.e. /3 = 0) and (PTjl becomes the standard 
dispersion equation of an extraordinary wave propagat- 
ing in an uniaxial dielectric mediimr [21]. Note that the 



6 



toroidal metamaterial does not interact with electromag- 
netic field, polarized perpendicular to n. 

Figure [5] illustrates the extent to which the effect of 
the mutual coupling between the toroidal metaparticles 
affects the dispersion properties of the wave propagating 
in a bulk sample of the metamaterial. As can be seen if 
the Q-factor ujo /j of an individual toroidal resonator is of 
the order or below 200, the mutual interaction between 
the metaparticles plays no role (Fig. [SKa)), although the 
toroidal medium in the latter case can be formally re- 
garded as "dense" (i.e. adjacent metapaticles have been 
brought in contact with each other and lot — ^to)- This 
is because, by its nature, toroidal interaction with exter- 
nal electromagnetic fields is much weaker than electric- 
and magnetic-dipole interaction with the same field. The 
effective electromotive forces driving the currents in the 
toroidal solenoinds are, therefore, relatively small and the 
presence of stronger losses further limits the magnitude 
of these currents. Thus, the mutual interaction between 
toroidal moments is a higher-order effect even for " dense" 
toroidal media. In contrast, at higher Q-factor values 
(Fig. Eljb)), the interaction between the metaparticles be- 
comes significant. Figure[S]compares the electromagnetic 
properties of a toroidal metamaterial (a), (c) to that of 
an array of split-ring resonators (b), (d). As can be seen, 
toroidal response shows a strong dependence on the back- 
ground dielectric permittivity e of the host material. In 
particular, doubling e effectively doubles the stop-band 
bandwidth of the toroidal medium (a) , whereas doubling 
e has little effect on the response bandwidth of an array of 
split-ring resonators (b). This feature suggests possible 
sensor applications, since e depends on parameters such 
as temperature, stress, external electric field, etc. Thus 
small relative changes of e may result in significant rel- 
ative changes of the transmission coefficient. The origin 
of this feature can be traced back to the magnetoelectric 
nature of the toroidal response. Since kx H = —ueeoE, 
Eq. (fT5)) suggests that the macroscopic magnetization 
M is effectively driven by the electric field E. Neglect- 
ing the interaction between the metaparticles, setting fc^ 
to zero and £p to 1 in (|2ip results in fc^ = ^j:^ with 

(the ratio be- 



The quantity 



I .£ 

I k ds 



tween the relative variation of the wavenumber k and the 
corresponding relative variation of the ambient permit- 
tivity e) can be used to estimate the "sensitivity" of the 
system. The result is 



e dk 
k de 



il + 2ReigT) + \9T\T^- 



The resonant structure of gri^^) suggests the possibility 
that near the resonance 1 -I- 2Re{gT) + Iffrl^ ~ and this 
results in a strong dependence of the power transmision 
coefficient |T(a;)p on e. It should be emphasized that the 
latter result is a feature that is unique for toroidal media 
and does not exist for ordinary dielectric/magnetic mate- 
rials composed of arrays of electric/magnetic dipoles. In- 
deed, since k^ = ^-e/i in an ordinary dielectric/magnetic 
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FIG. 6: Frequency dependence of tlie power transmission coef- 
ficient \T{ijj)^ for a metamaterial slab of tliickness /o = 10cm 
at normal incidence and iS||n. (a) Array of toroids with 
d = 5mm, R = d/3, ujt /2-k = 28.3GHz, ljo/2-k = 2.5GHz 
and ^/loo — 4.4 x 10"'^. (b) Array of split-ring resonators 

H •> •> — : — , where 



with an effective permeability ^J.{Lo) = 1 + jij'j. 

F — 0.1. The other parameter values are the same as in (a). 
In both (a) and (b) the full (red) line corresponds to e = 1 
and dashed (blue) line to e = 2. (c) Composite toroidal meta- 
material (see Fig. [Sjl . The wires are assumed lossless and the 
other parameter values are the same as in (a), (d) Composite 
metamaterial consisting of an array of lossless wires and an 
array of split-ring resonators with the same effective perme- 
ability as in (b) . In (c) and (d) the full (red) curve pertains to 
the case ujp — \/^u}o, dashed (blue) curve - to ujp = 1.5\/2ujo 
and "f/uo = 4.4 x 10"". 



material and, hence, gr — 0, the "sensitivity" |f ^| re- 
mains relatively low at all frequencies in full accordance 
with Figure inijb) and (d). The Q-factor ujq/j of a single 
toroidal metaparticle in Fig. ©(a) is 230, which is com- 
parable to the value of 50 measured by the authors of 
in the same frequency range, but for toroidal solenoids 
that are 10 times smaller than those considered here. At 
the same time the analysis indicates that higher Q-values 
should be achievable with larger solenoids [IjI- As Fig. 
©(c) shows the same strong dependence of the toroidal 
medium response bandwidth exists for a host medium 
with negative permittivity whereas no such feature is ob- 
served for arrays of split-ring resonators (d). In this case 
the stop bands from (a) and (b) are transformed into 
pass-bands by the presence of the negative permittivity 
material. Note, however, that much higher Q-factors are 
needed to produce a significant transmission in the com- 
posite toroidal medium. 

Backward waves exist at frequencies within the pass- 
band shown in Fig. ©(c). To appreciate this consider 
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(a) i (b) i 

FIG. 7: Refraction properties of a composite toroidal meta- 
material for a p-polarized incident wave. Poynting vector S 
and wavevector k of the incident ("i") and the transmitted 
("J") wave, (a) The optic axis n is parallel to the interface, 
(b) The optic axis is perpendicular to the interface. 

the Poynting vector S, in the absence of losses for fc_Ln 
and E\ \n 

S='^sil~ic.,/c.r)il+g{Lu)), (23) 

up' 4 

where o = ^ H — •! i 9/-. , ? Z*^ - i ^^ — 1\ and fc^ = 
^Y^. This shows that if 1 + giyi) > and oj < ujp 
the material is transparent and the wavevector and the 
Poynting vector are antiparallel. Analysis similar to that 



of [221 has been performed and the results are summa- 
rized in Fig. [7l The components of the Poynting vector 
and the wavevector that are perpendicular to the optic 
axis point in opposite directions - i.e. S±.k± < 0. De- 
pending on the orientation of the optic axis with respect 
to the interface this results in either positive or nega- 
tive refraction, shown on Fig. [Tl^a) and (b), respectively. 
Note, however, that despite the possibility of negative 
refraction (Fig. [Tfb)) the composite toroidal metamate- 
rial is not a left-handed metamaterial [2^ since a wave 
propagating in a left-handed metamaterial is always a 
backward wave - it satisfies the more restrictive condi- 
tion S.k < 0. 

In summary, it is shown that in a sharp contrast to 
materials that exist in nature, a new type of toroidal 
metamaterial shows a significant toroidal response. The 
magnetoelectric nature of the toroidal response results 
in strong dependence of the transmission properties on 
the dielectric permittivity of the host medium, which 
suggests possible applications. Negative refraction and 
backward waves are found in a composite toroidal meta- 
material. 
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